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Abstract 

We investigate exact results of isotropic turbulence in three-dimensions when the 
pressure gradient is negligible. We derive exact two-point correlation functions of 
density in three-dimensions and show that the density-density correlator behaves as 



|xi — X2I where as = 2 + -^^P. It is shown that, in three-dimensions, the energy 
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p) ^ spectrum E{k) in the inertial range scales with exponent 2 — ~ 1.5212. We also 

. discuss the time scale for which our exact results are valid for strong 3D-turbulence in 



the presence of the pressure. We confirm our predictions by using the recent results of 
numerical calculations and experiment. 

PACS numbers 47.27.AK, 47.27.Jv 
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1- Introduction 

Recently, tremendous activities have started on the non-perturbative understanding of 
turbulence [1-12]. A statistical theory of turbulence has been put forward by Kolmogorov 
[13], and further developed by others [15-17]. The approach is to model turbulence us- 
ing stochastic partial differential equations. The simplest approach to turbulence is the 
Kolmogorov's dimensional analysis, which leads to the celebrated k"^/^ law for the energy 
spectrum. This is obtained by decreeing that the energy spectrum depends neither on the 
wavenumber where most of the energy resides, nor on the wavenumber of viscous dissipation. 
Kolmogorov conjectured that the scaling exponents are universal, independent of the statis- 
tics of large-scale fluctuation and the mechanism of the viscous damping, when the Reynolds 
number is sufficiently large. In fact the idea of universality is based on the notion of the 
" inertial subrange" . By inertial subrange we mean that for very large values of the Reynolds 
number there is a wide separation between the scale energy input L and the typical viscous 
dissipation scale r] at which viscous friction become important and the energy is turned into 
heat. 

However, recently it has been found that there is relation between the probability distri- 
bution function (PDF) of velocity and those of the external force [18]. This observation has 
been confirmed by experiments [19], and numerical simulations [20]. 

In this direction, Polyakov [5] has recently offered a field theoretic method to derive the 
probability distribution or density of states in (l-l-l)-dimensions in the problem of randomly 
driven Burgers equation [21]. In one dimension, turbulence without pressure is described by 
Burgers equation (see also [14] concerning the relation between Burgers equation and KPZ- 
equation). In the limit of high Reynold's number, using the operator product expansion 
(OPE), Polyakov reduces the problem of computation of correlation functions in the inertial 
subrange, to the solution of a certain partial differential equation [22,23], see also [28], about 
generalization of Polyakov's approach, to find the probablity density and scaling exponent of 
the moments of " longitudinal" velocity difference in the three-dimensional strong turbulence. 

In this paper we consider three-dimensional isotropic turbulence without pressure, which 
is described by Navier-Stokes equations, when the pressure gradient is negligible. We derive 
the Polyakov's master equations in higher dimensions and solve it, in the three dimensions. 
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We derive the exact exponent of two-point density correlation functions and the energy 
spectrum exponent. We also discuss the time scale for which our exact results are valid for 
strong 3D-turbulence in the presence of the pressure. 

2- Turbulence Without Pressure in Three— Dimensions 

We consider the following quantity: 

eA = p(x, t) exp(A.u(x)) (1) 
where p and u are the density and the velocity satisfying the Navier-Stokes equations: 

ut + (u • V)u = i/V^u - ^ + f (x, t) (2) 

Pt + daipUa) = (3) 

where p and v are the pressure and viscosity, respectively. The stirring force f(x, is a 
Gaussian random force with the following correlation: 

< /^(x, i)/,(x', t) >= A;^,(x - x')5(t - t) (4) 

where iJ>,iy — Xi, X2, - ■ ■ , xn- 

We start with the situation when Vp ~ 0. This mode has a characteristic time of 7], ~ 
where L and Cg are the large scale of the system and the velocity of sound in the turbulent 
system, respectively. The existence of this time scale means that for the times t < Tp, we 
can ignore the pressure term in the Navier-Stokes equations and therefore our results are 
valid only for this time scale. 

To investigate the statistical description of eqs.(2) and (3), we consider the following 
two-point generating functional 

F2(Ai, Aa, xi, xa) =< p(xi)p(x2) exp(Ai.u(xi) + A2.u(x2)) > (5) 

We write the eq.(2) and (3) in two points Xi and X2 ior ui,U2, ■ ■ ■ ,un and p{x) and mul- 
tiply the equations in ^(xa), Ai^ip(xi)p(x2), • • •, Aia,^p(xi)p(x2) and p(xi), A2a;iP(xi)p(x2), 
• • •, and A2a;jvP(xi)p(x2), respectively. 

We add the equations and multiply the result to exp(Ai.u(xi) + A2.u(x2)) and make 
averaging with respect to external random force, therefore we find: 
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9tF2+ E dX~^'^^'^ ^ C,,, + D2 (6) 

{i=l,2}ij,=xi,---,XN ^'^ {i=l,2}/i=a;i,---,a;jv 

where Cj_^ and D2 are, 

E C'^-M^ E Ai,^Aj>/c^,.(xi -Xj)F2 (7) 

{i=l,2}/i=i;i,---,a;iv {i,j=l,2}ix,u=xi,---,XN 

and 

i:)2 =< i^p(xi)p(x2)[Ai • V^u(xi) + A2 • V^u(x2)] exp (Ai • u(xi) + A2 • u(x2)) > (8) 

where we have used the Novikov's theorem. D2 is known as the anomaly term [5] . 

Now we consider the anomaly term in the limit of small v or high Reynold's numbers. It 
is noted that this term can not be written in terms of F2. To find its structure we consider 
the symmetries of the basic equations. The basic equations are Galilean invariant and also 
are invariant under the rescaling of density as p — ap. In the other hand, final expression for 
D2 must contain the ultraviolet finite operators Vu, p and e'^'". The only finite combination 
satisfying the rescaling p — > ap is pc^'" (see [24] for more detail). Therefore D2 has the 
following form: 

D2 = aF2 (9) 

where a is generally a function of Ai and A2. 

Therefore in the steady state we find the following equation for F2: 

{i=l,2}ij,=xi,---,XN ^'/^ 

- E K,n>^j,uk^,u{xi - Xj)F2 = aF2. (10) 

{i,j=l,2}fi,u=xi,---,XM 

Also in this paper, we suppose that k^i, has the following form: 

v(xi - xj) = k{Q)[i - - (11) 

with k{0),L = 1. 

Now let us consider the eq.(lO) in three dimensions. We change the variables as: Xj. = 
Xi ± X2, A+ = Ai + A2 and A_ = ^^2^^ and and consider the spherical coordinates, so that 
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x_ : (r, 9, (p) and A_ : (p', 9', </?'). Direct calculation shows that: 
^ d d 

i=l,fj.=x,y,z fJ.=x,y,z '-'^-ti 



sin9cos9'cos{(p — </?') — cos9sin9' ^ ^ ^ sin9sin{(p — cp') 

p'sin9' 



+ ; drd0> H r^-^. drd^> 



sin9'cos9cos{p — p') — cos9'sin6 , cos9sin{p ~ 

+ o de H ^— — OeO^' 

r rpsin9' 

_ cos9'sin{p ~^') Q Q ^ ^^^iv d d . 

rp'sinO ^ ^ rp'sin9sin9' 

cos9cos9'cos(p — p') + sin9sin9' ^ „ 

+ , dedff, 

rp 

sin9'sin{(p - (p') ^ ^ 

Tsin9 '^'''^^ ^ ^ 



and 



5Z \,n^j,vk^i,{'^i Xj) 

{«J=l,2}At,i'=2^,y,2 

= [rV + 2(x_A_a. + y-\-y + z-\-,f] = rV"(l + 200527) (13) 

where cos^ = cos9cos9' + sin9sin9'cos{p — p'). 

Now using the eqs.(12,13) for isotropic turbulence we obtain: 

S(l-S2) 1 + ^2 

+ 

+ - rV'(l + '2s')]F^ = a{p')F^ (14) 

p r 

where C0S7 = s. The p' dependence of the a{p') anomaly must be chosen to conform 
the scaling and can be different depending on the scaling properties of the force correlation 
functions. In general, in the case of isotropic turbulence, stirring correlation function behaves 
as A;^,t/ ~ 1 — r^', where in our case wc have rj = 2. Therefore, a must depend on p' as follows 
a{p') — aop''^, where a — It is evident that for our case a is independent of p'. Let us 
consider the universal scaling invariant solution of eq.(33) in the following form: 

F2{p',r) = g{r)F{p'r^) g{r) = r'^^ (15) 

where 5 = and is the exponent of two-point correlation functions of density and also 
using the eq.(ll) we find 6 = 1. 



We substitute eq.(15) in eq.(14), and find the following relation for F{p'r): 

-cgi^a, + ^^d,dr + ^—^d,d' - rV'(l + 2s')]F{p'r) = aoF{p'r) (16) 

Since the two point density correlators exists, in the limit of p' 0, F{p'r) tends to a 
constant and thus we have to look for the solution of F among the family of positive, finite and 
normalizable solution of eq.(16). In the other hands, taking the Laplace transformation of the 
above equation, one can show that, to consider physical solution, so that < u(xi) — u(x2) >= 
0, we have to consider the case ao = 0. However for different types of correlation for the 
stirring force, e.g. k^^u ~ 1 — r'', with r/ 7^ 2, we have to include the oq [23]. 

Now, we propose the following Ansatz for F{p'r), with z — p'r: 

F{z, s) = e^'^(^) (17) 
Using the eq.(16), we find 7 = 3/2 and f{s) satisfy the following equations: 



Isfisf + 3f{s)f{s){l - s') + f{sy{-s + s') = (1 + 2s') 
-s{l - s')f{s)" + [(4 + as) - (2 + as)sV{s) + (f + |a3)^/(^) = 

also from eq.(18-a), one can derive the following initial conditions for f{s), 



(18) 



/(I) = -I /(I) = ^^i±^ (19) 

It is interesting to note that the equation for f{s) (i.e. eq.(18-a)) is the same as equation 
which is found in the instanton approach [8] . 

The function f{s) has the following expansion around s = 1: 

2 ^+\/n, , 5^33-61 , ,2 

- 75 ^ - ^ M - 2yn) - ■ ■ ■ 

Now using the boundary conditions on f{s) (i.e. eq(19)), and positivity of the probability 
distribution function we find: 

«3 = ^ 2.9574 (21) 
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Noting the fact that p' has dimension —1, we can find the following scaling relation for 
the density of the energy e{x\ 

e{ax) = a^e{x) (22) 

where A = 1 — and therefore we can determine the behavior of the energy spectrum 
exactly as: 

E{k) ~ k-f^, /3 = 2 - ^ ~ 1.52128 (23) 

This behaviour of energy spectrum is known as the non-Kolmogorov power laws which 
has been observed experimentally [25,27] and also in the numerical simulations [25,26]. 

The numerical calculations have been done in [25,26], where they have used the Wiener- 
Hermit expansion. They have shown that the energy spectrum behaves as E{k) ~ k^^'^'^^ for 
systems without boundaries (i.e. free turbulence) and also for a finite system this spectrum 
is not stable. In [25,26] it has been shown that in the inertial subrange for a finite system 
energy spectrum starts with slope —1.521, and after a moderate time which is less than 
the characteristic time ~ — ^ ( where L and Urms are the large scale of the system and 
rms value of initial velocity fluctuation, respectively) the equilibrium is attained and has 
transformed to —5/3. In other words the — 5/3's law is the stable algebraic spectrum for the 
Navier-Stokes equations after a time of order Tc. 

The experimental results (reported by Wissler [25,27]) shows that the non-Kolmogorov 
spectrum has been observed also experimentally only for moderate times less than Tf.. It is 
noted that in general for a turbulent flow Urms < Cs and therefore the pressure time scale 
has the property that Tp <Tc. 

Finally we can derive the PDF for the velocity difference and show that it tails as e~"" in 
the limit |w| — > +oo, which is in agreement with other approaches [18] for three-dimensional 
turbulence. 

Acknowledgements: We would like to thank B. Davoudi , J. Davoudi, R. Ejtehadi, M. 
Jafarizade, A. Langari M.R. Mohayaee and S. Rouhani for useful discussions and D.D. 
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